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THE  INSTABILITY  BOUNDARY  AND  GROWTH  RATE  CONTOURS  OF  THE 
GYROTRON  AMPLIFIER  WITH  SPREADS  IN  ENERGY  AND  MOMENTA 


1-3 

The  mechanism  of  the  electron  cyclotron  maser  instability  has  been 

4 

used  successfully  in  the  gyrotron  oscillator  which  provides  a  powerful  source 
of  millimeter  and  submillimeter  microwave  generation.  Many  reviews  on  the 
developments  of  theories  and  experiments  of  the  gyrotron  already  exist  in 
the  literature.^  ^ 

Recently  the  bandwidth  and  linear  growth  rate  of  a  gyrotron  amplifier  has 
been  investigated  in  theories^  ^  and  experiments. The  measured  linear 
gain  of  a  gyrotron  amplifier  lies  in  the  predicted  range  of  the  growth  rate. 

Its  poor  efficiency^  however  was  attributed  to  the  velocity  spread  effect, 
i.e.,  the  spreads  in  momentum  and  energy.  It  is  our  objective  to  investigate 
the  effects  of  the  spreads  on  the  linear  growth  rates  over  a  wide  range  of 
the  mode  spectra. 

We  consider  a  gyrotron  amplifier  which  consists  of  the  metallic  cylindrical 
waveguide  with  radius  R,  and  the  annular  tenuous  electron  beam  under  the  uniform 
axial  magnetic  field  Bq  inside  the  waveguide.  We  choose  the  equilibrium  distri¬ 
bution  function  of  beam  electrons  fQ,  given  in  terms  of  the  total  relativistic 
mass  factor  y,  the  axial  linear  momentum  pz,  and  the  azimuthal  canonical 
angular  momemtun  P. ,  to  be  a  rectangular  type  function  in  y,  p  and  Pa  with 

the  mean  values  given  by  y,  pz  and  PQ,  and  their  spreads  by 

Note:  Manuacript  submitted  September  11,  1979. 


1 


mode  is 


ye^,  p2ez  and  P0e0(Fig.  1).  The  linear  dispersion  relation  for  a  TEing 
then  obtained  for  the  eigenf requency  o>  at  the  azimuthal  mode  number  £■ ,  the 
radial  mode  number  n,  and  the  magnetic  harmonic  number  s,  as  a  function  of  the 
axial  wave  number  k. 

In  order  to  be  managable  we  assume  that  the  spreads  are  not  large  such 
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that  only  their  effects  on  the  resonance  terms  are  kept,  and  further  assume 
that  the  growth  rate  is  small,  i.e.,  |u>J  <<a>c/y,  where  <*>c  =  I  e|  BQ/mc.  Then 

the  dispersion  equation  is  reduced  to  the  equation  * 

A+a+  nbF  (A,  e  .  e  )  -  0,  (1) 

I  Z 

where  A =u-fts  represents  the  eigenf requency  shift  from  the  Doppler-shifted 

cyclotron  frequency  of  beam  electrons  ft  ,  defined  by  ft  =kc&  +su>  /y  ; 

s  s  z  c 

2  2 

a  H  (ftg  )/2ftg  refers  to  the  degree  of  the  synchronization  between  the 

2  2  2  2 

beam  mode  ft_  and  the  waveguide  mode  u.  defined  byu„  *c  k  +w  with  u 

8  in  1  in  co  co 

2  2  2 

as  the  cutoff  frequency;  b  =(ft  -c  k  )/2ft  measures  the  available  energy;  and 

s  s 
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n  is  the  beam-waveguide  coupling  coefficient.  Here  cBz  is  the  average 

axial  velocity  of  electrons.  The  small  growth  rate  assumption  (  |a|  <<  wc/y) 

_2 

favors  the  instability  term  (proportional  to  A  )  over  the  stabilizing  terms 
(proportional  to  A  ^  or  A°  )  in  the  dispersion  equation.  Consequently 
equation  (1)  slightly  overestimates  the  actual  growth  rate.  The  information 
on  the  influence  of  the  spreads  are  contained  in  function  F,  defined  by 

*Equatlon  (1)  vaa  first  obtainad  by  Uha  and  Davidson  (l.a.,  s««  Eq.  29  of  tha  paper;  Phys. 

Fluids  22,  September  1979)  for  the  TEon  node.  Here  we  generalised  the  case  to  the  azlomth- 
ally  asymmetric  TE.  mode.  The  analysis  Is  carried  out  In  this  report  such  as  to  keep  the 
combined  effect  of  both  spreads  maintained  for  the  arbitrary  axial  wave  number  k. 


2 


The  function  F  is  interchangeable  with 


respect  to  e.  and  e-,  and  an  even  function  of  e.  and  e_ ;  F(A,e 


recovers  the  dispersion  relation  without  spread 


and  e  are  weighted  by  |  f2  |  and  |  kcB  |,  we  find  that  the  influence  of  the 

2  S  Z 

energy  spread  is  dominant  over  that  of  the  axial  momentum  for  small  |k|,  while 


for  large  |k|  the  reverse  is  true.  In  usual  gyrotron  devices,  the  wavenumber  k 


is  commonly  chosen  such  that  and  e,  are  in  the  same  order  of  magnitude 


The  growth  rate,  viz.  the  imaginary  part  of  A,  depends  on  several  experi 
mental  parameters:  V  (beam  voltage),  I  (beam  current),  g  /g  (ratio  of  the 


velocity  in  perpendicular  to  axial  direction),  r  (beam  center  location) 


n,  R  and  k  in  addition  to  e  and  e  .  By  introducing  a  normalization 


equation 


the  dispersion  equation  (1)  becomes 


The  advantage  of  normalization  (3)  becomes  apparent  when  we  note  that  the  re 
duced  dispersion  equation  (4)  now  becomes  a  function  of  mere  dimensionless 


numbers 


and  that  the  experimental  parameters  are  absorbed 


in  normalizing  factor  (nb) 


•  pjM 


as  ^unct^ons  the  wave  number  k  are  Illustrated  in  Fig.  2.  It  shows  a 

TEnil  mode  with  V  ■  71.5  kV,  I  ■  9.5A,  6^/8  ■  1.5,  and  r  and  R  optimized 

to  give  highest  growth  rates.  ^  Note  that  a  £  0  for  and  the  upper 

limit  of  a  is  determined  by  the  wall  radius  R.  We  also  note  that  the  scaling 
1/3  -2 

factor  (nb)  is  in  the  order  of  10  and  remains  relatively  unchanged  over 

k  except  for  the  stable  regions  at  the  edges  in  Fig.  2.  The  weighting  factor 

for  e  is  much  larger  (''.AO)  than  that  for  e  C'.3).  Therefore,  in  order  to  have 
Y  2 

equivalent  effect,  only  a  fraction  of  the  momentum  spread  (e  <  0.1  e  )  is 

Y  2 

sufficient  for  the  energy  spread. 

Examination  of  the  dispersion  equation  (A)  yields  that  there  exists  an  upper 
limit  of  the  growth  rate,  given  by  |^|  <  min  {3*4/2,  ir/Ae^t^}.  The  absolute 

<\j  Ik  f\j  'Xj 

maximum  growth  rate  ■  3  /2  occurs  for  the  zero  temperatured  (e^  «  e2  «  0) 

•\j 

electron  beam  at  the  beam-waveguide  synchronous  (0  *  w.  )  mode  (a  *  0) . 

s  fen 

This  shows  that  the  rate  of  energy  extraction  from  beam  electrons  is  highest 
when  the  beam  and  the  waveguide  modes  are  synchronous.  The  possible  maximum 

<\,  r\, 

gain  reduces  as  the  spread  (e^  or  e2)  increases.  Equation  (A)  also  shows  that 

■V.  'Xi  *\j 

the  instability  region  is  bounded  in  a  a-space  by  a+  >a  >a_.  The  instability 

'Xj  %  'x, 

boundaries  a+  and  a_  are  given  in  Fig.  3.  The  lower  bound  a  is  a  decreasing 

function  of  £2(or  e^)  ,  while  the  upper  bound  a+  is  a  decreasing  function  of 
|e^  -  e2l  with  its  maximum  at  *  e2-  Since  the  upper  limit  of  a  is  experi¬ 

mentally  imposed  by  wall  radius  R  (see  Fig.  2)  unless  the  spreads  are  very 
large,  in  which  case  it  is  determined  by  a+  it  is  found  that  for  small  e  the 
unstable  range  in  k  is  determined  by  a_.  Therefore  we  conclude  that  the  unstable 
bandwidth  in  the  wave  number  k  becomes  broader  as  the  spreads  increase. 


Finally  we  proceed  to  find  Che  eigenvalue  X  by  numerically  solving  eq.  (4) 

f\j  f\j  % 

for  given  a,  and  £2*  The  results  are  shown  in  Fig.  4.  The  growth  rates 

f\j  *\j  *v> 

contours  are  drawn  in  e.-e-  space  for  four  representative  values  of  a  (1,0,-1 


The  instability  region  in  spread  space  becomes  broader  as  a  decreases  due  to 


0  due  to  a_  curve.  The  maximum  growth  rate  in 


0  in  general,  but  as  a  decreases,  the  point 


marked  by  a  filled  square,  moves  along  t 


This  phenomenon  is  readily  explained  when  we  note  that  at  out  of  synchronization 


mode  (a  *0)  the  spreads  provide  a  fraction  of  electrons  that  are  synchronous 


while  at  the  synchronous  mode  (a  ■  0)  the  spreads  reduce  the  number  of  synchro 


nous  electrons.  Gyrotron  amplifiers  are  usually  operated  near  a  »  0.  With 


parameters  given  in  Fig.  2,  the  minimum  spreads  needed  to  stabilize  the  insta 


bility  are  found  to  be  3.4%  and  49.2%  for  the  energy  and  the  momentum  at  the 


synchronous  case.  Near  the  beam-waveguide  synchronous  mode  (i.e 


while  the  growth 


the  instability  is  stabilized  with  the  spreads  for  low  e 


rate  is  maximum  at  e 


synchronization  mode  (i.e 


>>  0).  The  opposite  behavior  of  the 


growth  rate  on  the  spreads  at  high  e  may  be  significant  if  the  spreads  increase 


during  the  interaction  stage 


We  note  that  at  the  synchronous  mode  the  growth  rate  becomes  highest  along  the 


0).  At  the  synchronous  condition  where 


the  conducting  wall  radius  R  is  adjusted  to  have  the  beam  mode  ft  curve  graze 


corresponds  that  6^.  e 


the  waveguide  u.  curve,  the  line  e 


where  e  refers  to  the  spread  in  the  transverse  momentum.  This  means  that 


I 

the  thermal  temperature  of  the  beam  electrons  Is  Isotropic.  In  other  words, 
at  the  synchronous  wave  number  the  amplifier  becomes  most  efficient  when  the 
beam  distribution  is  isotropic. 

It  is  worthwhile  to  emphasize  that  the  reduced  dispersion  equation  (4) 
coupled  with  the  normalization  equation  (3)  describes  the  entire  spectrum  of 
mode  s,  £,  n,  and  k.  Once  the  mode  is  identified,  eq.  (3)  determines  the 

'X,  'Xj  *\j  y 

numerical  values  of  a,  and  tj  to  find  the  growth  rate  via  eq.  (4), 
which  in  turn  is  transformed  to  the  actual  through  eq.  (3).  Therefore, 
the  mode  spectral  analysis  is  reduced  to  determination  of  the  scaling  factor 
eq.  (3). 

The  authors  would  like  to  acknowledge  the  benefit  of  the  useful  discussions 
with  H.  Uhm. 
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Fig.  2  -  Doppler  shifted  beam  mode  (  ftg)  and  the  waveguide  mode 
vs.  the  axial  wavenumber  k,  and  the  corresponding  normal¬ 
ization  factor(nb)l/3^  C^e  norma^^ze^  a  vs*  ^ 

For  actual  parameters  see  text. 
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Fig.  3  -  Instability  boundary.  The  mode  is  unstable 
for  the  range  k  <  k  <  k+. 
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